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^ : Abstract 

^^ . The problem of restoring Proissart bound to the BFKL-Pomeron is studied in 

an extended leading-log approximation of QCD. We consider parton-parton 

> 

0^ . scattering amplitude and show that the sum of all Feynman-diagram contribu- 

(N 

O 

Q^ . by the phase shift, and subleading-logs of all orders needed to restore the 

r~| ' Froissart bound are automatically provided. The main technical difficulty is 

JL , to find a way to extract these subleading contributions without having to com- 

,^ . pute each Feynman diagram beyond the leading order. We solve that problem 



tions can be written in an eikonal form. In this form dynamics is determined 



by using nonabelian cut diagrams introduced elsewhere. They can be con- 
sidered as colour filters used to isolate the multi-Reggeon contributions that 
supply these subleading-log terms. Illustration of the formalism is given for 
amplitudes and phase shifts up to three loops. For diffractive scattering, only 
phase shifts governed by one and two Reggeon exchanges are needed. They 
can be computed from the leading-log- Reggeon and the BFKL-Pomeron am- 
plitudes. In applications, we argue that the dependence of the energy-growth 
exponent on virtuality Q^ for 7*P total cross section observed at HERA can 
be interpreted as the first sign of a slowdown of energy growth towards sat- 
isfying the Froissart bound. An attempt to understand these exponents with 
the present formalism is discussed. 



I. INTRODUCTION 

Large rapidity-gap events observed recently at HERA |l],0 provide additional impetus to 
QCD calculation of diffractive scattering and total cross section. At virtuality Q^ ^ ^'qcdi 
the QCD fine structure constant as = g"^ /^t^ is small, near-forward parton-parton scattering 
is in principle calculable in perturbation theory, and total cross section can be obtained via 
the optical theorem. To be sure, such calculation is very complicated at high cm. energy 
^/s, as the effective expansion parameter for the problem is a^lns and not just a^. The 
former could be large at high s even though the latter may be small. As a result, multi-loop 
diagrams must be included. 

To make multi-loop calculations feasible, leading-log approximation (LLA) is usually 
employed. This consists of keeping only the highest power of In s at each order. Equivalently, 
if the scattering amplitude (divided by s) is considered as a function of the two variables a^ 
and Us In s, then only the lowest power of as is kept. When the resulting amplitude is of the 
form a"sF(aslns) for some function F, we will abbreviate it simply as ~a". 

The LLA result for high-energy near-forward scattering is known, though not completely. 
The dominant parton-parton scattering amplitude occurs with the exchange of a colour-octet 
object called the Reggeon (or the Reggeized gluon), and is ~a;s 0. The Reggeon amplitude 
is actually proportional to {ags/ /S?')s^^^\ with A^ = —t the square of the momentum 
transfer, and -R(A) a known function of A proportional to a^. Note that s^^^^ is of the form 
^(aslns), so the Reggeon amplitude is indeed ~as, as claimed. Although the Reggeon 
amplitude is the dominant result for summing all diagrams in LLA, it can be obtained just 
from the t-channel ladder diagrams and their gauge partners, because all the other diagrams 
are subdominant. We may therefore think of the Reggeon as a composite object made up 
of gluons, bundled up together roughly in a ladder configuration. 

Diffractive scattering occurs via a colour-singlet exchange. Two or more Reggeons are 
required to form a colour-singlet object in the t channel. Since each Reggeon carries a small 
factor ^tts, the dominant singlet amplitude comes from two interacting Reggeons, giving 



a magnitude ~af. This is the BFKL Pomeron [g^Q. Unfortunately this amphtude is not 
known too precisely, though one does know that in the forward direction and extremely high 
energies it behaves like s'^ , with J > 1 a known number. This gives rise to a total cross 
section with energy variation s'^"^, which violates the Froissart bound (Ins)^. 

Present data from HERA [|T],3 on 7*P total cross section are consistent with a power 
growth in s, though with an observed exponent smaller than J — 1 computed from the 
BFKL Pomeron. To satisfy unitarity and obey the Froissart bound, the exponent must 
eventually decrease to zero at very high energies, and that can be achieved theoretically 
only by including subleading-log contributions that have been hitherto neglected. The first 
subleading correction to the BFKL Pomeron have been computed [§], with the encouraging 
result of a smaller J, but also with a pathological behaviour that has not yet been fully 
resolved. In any case the Froissart bound is still violated. 

From the s-channel unitarity relation 2Im(T/-j) = I]„Tj„T„i, one expects all multiple- 
Reggeon exchanges to have to be included before we can restore full unitarity and the 
Froissart bound. As each additional Reggeon brings in an extra factor of ~as, they provide 
the required subleading-log contributions to all orders. The purpose of this paper is to discuss 
a formalism whereby this scenario can be implemented in pQCD. For other approaches see 
for example P|. 

To clarify what we have in mind, let us first pretend the partons to be scattered through 
a two-body instantaneous potential V{x). Born amplitude alone violates unitarity, but that 
can be fixed up by including all higher-order corrections. For small V, these higher-order 
terms are the analogs of the subleading contributions discussed above. These amplitudes 
can be summed up at high energy into an eikonal formula (see eq. (|l|)) where unitarity is 
restored. Dynamics is now specified by the phase shift 6{b), which is a function of the impact 
parameter b, and is linear in V. Higher-order amplitudes can be recovered by expanding the 
exponential into powers of the phase shift. 

A very similar scenario emerged if the two partons interacted through a multiple exchange 
of photons. In that case the interaction is no longer instantaneous, and cross diagrams are 



present. In terms of the unitarity relation, this means photons are present in the intermediate 
states and the inelastic channels in which they are produced must be included in the unitarity 
sum. However, the eikonal form can still be established with the help of the 'eikonal formula' 
P), and the phase shift is still given by the Born approximation amplitude (the Coulomb 
phase shift) g|Tg. 

There is nothing new about using the eikonal form to implement unitarity. In fact, a 
very large amount of phenomenology of two-body scattering has been carried out in that 
framework. The modern challenge is whether, and if so how, pQCD amplitudes can be uni- 
tarized that way. For on-shell near-forward hadronic scatterings, confinement is presumably 
important and we may not even be able to use pQCD. However, 7*P total cross section at 
large virtuality is expected to be calculable in pQCD, and the BFKL Pomeron is an attempt 
to do so. Unfortunately it violates the Froissart bound, so the implementation of unitarity 
in pQCD remains unsolved. This is the area where we hope to make some progress on. 
The main difficulty comes from the fact that terms with all powers of a^ (with fixed a^ In s) 
are needed to restore unitarity, but these are small terms extremely difficult to obtain from 
Feynman diagrams in the calculable regime a^ -C 1. 

The following physical picture [|TT| may help to visualize how the Froissart bound is 
restored. The rise of total cross section |13| predicted by the BFKL Pomeron @] may be 



attributed to an increased production of gluon jets at higher energies. When the energy 
gets really high, there are so many gluons around that they tend to overlap one another. 
When that happens, coherent effect becomes important and a destructive interference sets 
in to reduce the power growth, to a rate that eventually satisfies the Froissart bound (Ins)^. 
This mechanism suggests that whatever formalism we use to restore the Froissart bound, 
interference effects ought to be a central part of it. The 'nonabelian cut diagrams' we propose 
to use to solve this problem fits into this category, because they can be viewed as a way 
to organize the summation of Feynman diagrams to heighten the Bose-Einstein interference 
effects of identical gluons [|I^ . 



Returning to the formal mechanism to unitarize QCD, the first thought would be to 



try to imitate potential scattering and QED by including multiple exchanges of Reggeons 
and/or Pomerons. This is indeed the general idea but there are non-trivial problems to 
be solved with this approach. Unlike the photon in QED or the potential, Reggeon and 
Pomeron are composite objects, themselves made up of gluons and possibly quark pairs. It 
is therefore not clear whether we are allowed to exchange them as if they were elementary. 
For example, should we include diagrams where Reggeons are crossed? Or equivalently, can 
Reggeons be produced from the partons so that inelastic channels involving their production 
must be included in the unitarity sum like QED? Can these composite objects overlap and 
merge? If an effective theory of Reggeon/Pomeron equivalent to QCD existed without 
gluons, then presumably we would be able to answer all these questions once its precise 
dynamics is known. However, since gluons can be produced off Reggeons and Pomerons, 
they must be included in any effective theory [jl2[, in which case it is hard to exclude them 



from being exchanged as well. If we have to exchange gluons anyhow, we might as well 
go back to the original QCD theory which exchanges nothing but gluons and quarks, and 
whose dynamics is precisely known. At least we can avoid potential double counting if 
we do it that way. Nevertheless, gluon self-interaction and non-commutativity of colour 
matrices have no analog in the potential or the QED problem, so it is not clear how the 
very complicated QCD Feynman diagrams can be handled. It is even less certain that the 
results can be mimicked by the exchange of Reggeons and/or Pomerons which we hope to 
see. Worst of all, Feynman diagrams can hardly be calculated beyond LLA, but we need 
subleading contributions to build up the eikonal form and unitarity. How can we possibly 
get them without going through the impossible task of computing Feynman diagrams to 
highly subleading orders? 

Fortunately these questions can be answered and difficulties overcome, if we use the 
nonabelian cut diagrams introduced elsewhere [p!5|-p!8[] instead of the conventional Feynman 
diagrams. Reggeons and Pomerons emerge naturally, and phase shifts can be calculated 
using only the leading-log approximation. Due to the possibility of gluon production from 
Reggeons, phase shift is generally no longer linear in the (Reggeon) exchange, as was the 



case in potential scattering and QED. 

Nonabelian cut diagrams will be reviewed in Sees. 3 and 4. They are Feynman diagrams 
with slightly different 'Feynman rules'. A permuted sum of Feynman diagrams can be 
shown to be equal to a permuted sum of (nonabelian) cut diagrams, which is why we may 
calculate amplitudes using either of them. The main advantage of cut diagram is that 
identification with Reggeon becomes natural. In fact, each cut amplitude is given by a 
product of Reggeon (fragment) amplitudes. This factorization allows the eikonal form to be 
built up, and subleading terms to be computed just in the leading-log approximations. 

We have used the word Reggeon fragment amplitudes to mean finite-order amplitudes 
whose sum builds up the full Reggeon amplitude. Hence the word 'fragment'. There may be 
many distinct fragments even at a given order. For brevity, the word 'fragment' will often 
be dropped. 

In Sec. 2, the impact-parameter representation of a high-energy two-body amplitude is 
reviewed. A physical argument is given to show how the Froissart bound is restored when 
the bound-violating Born approximation is iterated and summed up into an eikonal form. 
In Sec. 3, nonabelian factorization formula |]T3|,|TB[ and nonabelian cut diagrams fT^'H ^°^ 



tree amplitudes are reviewed; nonabelian cut diagrams are simply Feynman diagrams with 
the factorization formula built in. It is important to note that this factorization has nothing 
to do with the usual factorization of hard physics from soft physics 0. The present one 
goes along the s-channel, whereas the usual one goes along the t-channel. In Sec. 4, the 
technique is applied to two-body scattering amplitudes. The advantage of using cut diagram 
over Feynman diagram to compute a permuted sum is discussed. The connection with 
Reggeons is identified. The question of factorizability of these amplitudes into irreducible 
parts is discussed in Sec. 5. This factorization allows the eikonal form to be built up and 
phase shifts computed from the irreducible amplitudes. Diffractive scattering is considered 
in Sec. 6. In that case phase shift may be restricted to those obtained from one and two 
irreducible Reggeon exchanges, and can be computed from the LLA-Reggeon and BFKL- 
Pomeron amplitudes. By using phase shifts in the eikonal form unitarization of the diffractive 
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amplitude is achieved. To illustrate the details of previous sections, quark-quark scattering 
amplitudes and phase shifts are given to the three-loop order in Sec. 7. Sec. 8 is devoted to 
the observed 7*P total cross sections. We argue that the dependence of the energy-growth 
exponent on virtuality Q^, as well as the size of the BFKL exponent compared to the 
observed ones, can be taken as indications that the softening of energy growth towards the 
Froissart-bound limit is already happening. Correct energy dependence at various virtuality 
Q^ can be reproduced with the formulae of Sec. 7. Finally, Appendix A contains some 
colour-algebra calculations. 

II. IMPACT-PARAMETER REPRESENTATION, PHASE SHIFTS, 
AND THE FROISSART BOUND 

Let A(s, A) be a parton-parton scattering amplitude at centre-of-mass energy y^ and 
momentum transfer A = | A|. At high energies A is transverse, and conjugate to the impact- 
parameter b. Bold letters like A and b are used to describe vectors in the transverse plane. 

The two-dimensional Fourier transform of A(s, A) defines the impact-parameter ampli- 
tude .4.(s, 6) by the eikonal formula 

A(s, A) = lis f (fhe'^^-^Ai-s, b) = lis f (fhe'"^-^ (l - e^*^^"'^)) . (1) 

For large s ~ Ak"^, angular momentum is given by / = kb, and 6{s, b) is just the phase shift 
at that angular momentum. For QCD parton scatterings, the initial and final partons may 
contain different colours, so the amplitudes A(s, A) and ^(s, 6), as well as the phase shift 
(5(s, 6), should be treated as colour matrices. Only the diagonal matrix elements arising from 
a colour-singlet exchange are truly elastic amplitudes. We shall denote them by Ai(s, A) 
and Ai{s^ b). In terms of them, the total cross section (Jt{s) is given by the optical theorem 
to be 

aT(s) = -Im[Ai(s,0)]. (2) 

s 



No approximation has been made to arrive at (P, so unitarity is exact and the Froissart 
bound is satisfied. If 6{s, b) is small, the exponential may be expanded and only the term 
linear in 6{s, b) kept. This gives the 'Born amplitude' and allows the phase shift to be 
thought of as the 'potential'. Thus if the interaction causing the scattering has a range 
/i~^, then we expect S{s,b) ~aexp(— /i6). At large impact parameters, a may be taken to 
be approximately independent of b. In the region ^b ^ 1, the phase shift is small and the 
impact amplitude A{s, b) in (|l]) ceases to contribute. The effective radius -R(s) of interaction 
may therefore be estimated from the condition (5(s, R) ~ 1. If a is an increasing function of 
s, then R{s) and hence ar ~ 7TR{sy also increase with s. In particular, if a ~ s"'"^ has a 
power growth, caused for example by the exchange of a spin J > 1 Pomeron or elementary 
particle, then R{s) ~ [(J — l)//i] Ins, and ax ~ (InsY, which is the Froissart bound. This 
conclusion is general and is qualitatively independent of the magnitude of J and the range 
fi~^. Unless a increases faster than a power of s, the Froissart bound is guaranteed when 
total cross-sections are calculated from the phase shift. On the other hand, if we used the 
'Born approximation' throughout, then the total cross section would have grown like s'^~^ 
all the way, and the Froissart bound would be violated. This is essentially what happens to 
the BFKL Pomeron [Q . The cure, very roughly speaking, is to use it as a phase shift rather 
than an amplitude. This is the general idea but details are more complicated. They will be 
discussed in Sec. 6. 

To make use of this unitarization mechanism we must find a way to calculate the phase 
shift. This is not simple for several reasons. According to (|ID, even when S{s, b) is computed 
just to the lowest order, the resulting A{s, A) contains terms of all orders. This suggests 
that a proper understanding of phase shifts cannot be obtained until we know how to sum 
an infinite number of Feynman diagrams. To put it differently, we have to learn how to deal 
with 2i6{s, b) = ln(l — A{s, b)), which consists of an infinite sum of all powers of the matrix 
amplitude A{s, b). 

Moreover, according to (|I|), to be successful a certain product structure must emerge 
out of the sum. If the phase shift is expanded in powers of the coupling constant g^, 
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S{s,b) ~ Z^m fl'^"^^''™'^ then the (2?7,)th order contribution to A{s,b) is given by a sum of 
products of the phase shifts 5^"^^\ with '}2'mi = n. Individual Feynman diagrams certainly 
do not factorize in this manner, and it is not immediately clear why sums of Feynman 
diagrams have this structure either. But unless we can get the sum into this factorized form 
there would seem to be no simple way to extract the phase shift from Feynman diagrams. 

In the next two sections we shall lay down the foundation which enables us to build up 
such a factorized form. 

III. EIKONAL APPROXIMATION AND FACTORIZATION OF NONABELIAN 

TREE AMPLITUDES 

The decomposition of sums of Feynman diagrams into sums of products of Reggeon 
amplitudes, or of phase shifts, is based on a nonabelian factorization formula [^,|l^. This 
formula for tree amplitudes is the nonabelian generalization of the eikonal formula f^ used 
in high-energy scattering in QED. It can be conveniently embedded into Feynman diagrams 
to turn them into nonabelian cut diagrams ||T7|J18[| . These items will be reviewed in the 
present section. 

The formula deals with a sum of tree amplitudes like Fig. 1, in which a fast particle of 
final momentum p"^ emerges from an initial particle with momentum p^ after the emission of 
n gluons of momenta k^. Spins and vertex factors are ignored in this section but they will be 
incorporated later. The gluons concerned need not be on shell. This allows the tree diagram 
to be a part of a much larger Feynman diagram, so that formulae developed for trees here 
are also useful for more complicated diagrams later. We shall assume the energy p^ ~ p' 
of the fast particle to be much larger than its mass m, its transverse momentum A = |A|, 
the square-root virtuality Q in case it is present, and all the components /cf of every gluon 
momentum; in short, larger than any other energy scale involved. It is convenient to use a 
somewhat unconventional system in which the final particle moves along the z-axis. In the 
light-cone coordinates defined by A^ = {A^ ± A^), the final four- momentum of the energetic 
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particle can then be written as p'^ = {p'^,p~ , p) = (^/i, 0, 0), where terms of o(l) have been 
dropped. The initial momentum p'^ of the fast particle then carries a transverse component 
A = I]r=i ^i- ■'■^ ^^^ eikonal approximation outlined above, 

{p' + Ky - m'^ + ie 2p' ■ K + K^ + ie ~ 2p' ■ K + ie ~ y/s K' + ie ^^ 

can be made on all propagators of the fast particle, provided K is the sum of any number of 
the /cj's. The factor s~2 in (^ is irrelevant for the rest of this section so it will be dropped 
and the propagators taken simply to be {K~ + ie)~^. 

The tree amplitude in Fig. 1 is then given by the product of a momentum factor a[cr] 
and a colour factor t[cr], where 

(n \ n—1 -| 

t[a]=tit2t3---tn. (4) 

The (5-function in a[a] is there to ensure the initial-state momentum p = p' + J27=i k to be on 
shell, for in the eikonal approximation its square is m'^ + \^J2'i=i ^i ■ Since this (5-function 
is not explicitly contained in the T-matrix, it should be removed at the end, but for the sake 
of a simple statement in the factorization formula (|^) below it is convenient to include it in 
a[a] for now. 

The colour factor t[cr] is given by a product of SU{Nc) colour matrices ta in the repre- 
sentation appropriate to the fast particle. 

Since gluons obey Bose-Einstein statistics, we must sum over all their permutations to 
obtain the complete tree amplitude. Let [a] = [aia2- ■ ■ o"„] indicate the ordering of gluons 
along the fast particle from left to right. This symbol will also be used to denote the 
corresponding tree diagram. In this notation Fig. 1 corresponds to [a] = [a] = [123 •■■n]. 
The tree amplitude for the diagram [a] will be denoted by a[(T]t[(T]; they are given by (H) 
with appropriate permutation of the gluon momenta. The complete tree amplitude is given 
by the sum of individual tree amplitudes over the n! permutations [a] of the permutation 
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group Sn- The factorization formula [0,0 states that this permuted sum can be replaced 
by a similar sum of the nonabelian cut amplitudes: 

E «MiM = E «N^[<]- (5) 

[o-]gS„ [CT]e5„ 

Just like the amplitudes on the left-hand side which can be obtained from Feynman diagrams, 
the (nonabelian) cut amplitudes a[crc]t[cr^] on the right-hand side can be obtained from 
nonabelian cut diagrams ||17| , [T8[| . The (nonabelian) cut diagrams are Feynman diagrams 
with cuts inserted at the appropriate propagators. The 'Feynman rules' for cut diagrams 
are the usual ones except at and around the cut propagators, as we shall explain. The 
position of the cuts depends on [a], and is given by the following rule. A cut is placed just 
to the right of gluon ctj if and only if ctj < aj for all j > i. We shall label a cut by a 
vertical bar, either on the diagram itself or in the corresponding permutation symbol [a]. 
The resulting nonabehan cut diagram will be denoted by [(tJ- For example, if [a] = [32145], 
then [aj = [321|4|5]. If [a] = [12354], then [aj = [1|2|3|54]. 

The momentum factor of a cut amplitude, a[ac\, is obtained from the momentum factor 
a[cr] of the Feynman amplitude by replacing the Feynman propagator {K~ +ie)~^ of every cut 
line by the Cutkosky propagator —2TTi6{K^). There is then a superficial similarity between a 
Cutkosky cut diagram and a nonabelian cut diagram. However, they are completely different 
and the cut diagrams referred to in this paper are exclusively nonabelian cut diagrams. 

From (H), it follows that a[ac\ is factorized into products of a's separated by cuts, which 
is why the formula is called the factorization formula. For example, if [ad = [321|4|5], then 
a[ac\ = a[321]a[4]a[5]. 

We mentioned below eq. (H) that the ^-function appearing in the amplitude must be 
removed at the end. This can be carried out by putting in Cutkosky propagators only where 
a vertical bar occurs. The overall 5-fucntion for the sum of the ' — ' components of momenta 
then disappears because no cut is ever put after the last entry cr„ of [a] = [o"icr2 ■ • ■ (Jn]. 

The complementary cut diagram (or c-cut diagram for short) [cr^] of a cut diagram [ad is 
one in which every cut line in [ad becomes uncut, and vice versa. If no cuts appear in [a'd, 
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the colour factor t[a'^ is simply the product of colour matrices, as in a Feynman diagram. 
When cuts are present, the product straddling an isolated cut is replaced by its commutator. 
If consecutive cuts occur, then the product is replaced by nested multiple commutators. For 
example, if [a] = [32145], then [aj = [321|4|5], [<] = [3|2|145], and t[<] = [tg, [t2,ti]]Ut5. 

For n = 2, eq. (|) reads a[12]t[12] + a[21]t[21] = a[l|2]t[12] + a[21]t[2|l] = a[l]a[2]tit2 + 
a[21][t2, ti]. This can easily be checked by direct calculation. Explicit check for n = 3 is also 
possible |jl5|, but for larger n a direct verification becomes very complicated. 



The factorization formula (|]) is combinatorial in nature, and is true whatever the matri- 
ces ti are. In particular, if all the matrices t, commute, as in QED, then the only surviving 
term on the right-hand side of (^ is the one where no commutator occurs, which means 
that [a'^] contains no cuts, and [ad has a cut at each propagator. This implies factorization 
of the tree amplitude (H) into a single product Y[7=i '^[^l; which is the usual eikonal formula 



PJlOll used in QED scatterings. It is this factorization that allows exponentiation to occur 
and the phase shift to be computed. 

In QCD, colour matrices do not commute, 

[ta,tb] = if abate, (6) 

but their commutators generate other colour matrices. This allows the nested multiple 
commutator of colour matrices to be interpreted as a source for a new adjoint-colour ohiect, 
which we call a Reggeon fragment. For parton-parton scattering in the weak coupling limit, 
they turn out to be the constituents of the Reggeon as we know it from the Regge-pole theory, 
hence the name. The word 'fragment' is there to clarify that this object is just part of the 
Reggeon, not the whole thing, but for brevity this word is often dropped. This algebraic 
characterization of a Reggeon is valid even in the strong coupling limit, or physical situations 
other than two-parton scatterings. Unlike a gluon which is a point-like particle, a Reggeon 
is an extended object made up of a bundle of gluons, each interacting with the energetic 
particle at a different point. Later we may also introduce interactions between gluons to tie 
them together. Note that there are at least as many distinct Reggeon fragments as there 
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are nested multiple commutators. With this interpretation, the cut amplitude a[crc]i^[o"c] is 
nothing but a product of M such Reggeon amplitudes, with M — 1 being the number of cuts 
in [crj- 

Two final remarks. First, the rule of inserting cuts explained above depends on how the 
gluons are labelled, though at the end of the calculation it should clearly not matter how it 
is done. We may adopt the labelling which is most convenient for our purpose. Secondly, we 
have assumed the fast particle to carry a large '+' component of the light-cone momentum, 
but we may equally well construct cut diagrams if it had a large ' — ' component instead. For 
parton-parton scattering, both are required. 

IV. PARTON-PARTON SCATTERING AMPLITUDES 

We will study two-body amplitudes in this section using cut diagrams. Cut diagrams 
are easier to compute than the corresponding Feynman diagrams. In s cancellations occuring 
in permuted sum of Feynman diagrams do not happen in this case, and they are directly 
related to the Reggeon (fragment) amplitudes. These are some of the advantages of using 
cut diagrams. 

Fig. 2 depicts a QCD Feynman diagram for the scattering of two energetic partons. 
The upper one carries an incoming momentum p'^ = {y/s, 0, 0) in light-cone coordinates, 
and the lower one carries an incoming momentum p2^ = (0, y/s, 0). Suppose there are rii 
gluons attached to the upper parton and n2 to the lower, then the ni\n2\ diagrams obtained 
by permuting the position of attachment forms a permuted set of diagrams. The sum of 
amplitudes for diagrams in the permuted set will be called a premuted sum of amplitudes. 
If we replace the Feynman tree attached to each energetic parton by the corresponding 
cut tree, as in Sec. 3, then we have a cut diagram. According to (^), a permuted sum of 
Feynman amplitudes is equal to a permuted sum of cut amplitudes, which is why we may 
use cut diagrams instead of Feynman diagrams for computation. However, when both parton 
trees are permuted we may overcount, in which case the final result has to be divided by a 
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symmetry factor S. This for example will happen to a permuted sum of s-channel ladder 
diagrams. If n gluons are exchanged between the two energetic partons, then there are 
n! distinct diagrams by permutation. The permuted set however contains (n!)^ diagrams, 
because each distinct diagram appears n! times in the set. Hence a symmetry factor S = n\ 
is required for this case. 

The scattering amplitude depends on the centre-of-mass energy y/s and the coupling 
constant g, as well as the momentum transfer A = |A|, and the virtuality Q^ if it is 
involved in deep inelastic scatterings. Let us concentrate on its dependence on g and s for 
the moment. The Born amplitude is g'^s. An £-loop Feynman diagram is proportional to 
g2{i+i)^ but each loop may also (though may not) produce a factor of In s upon integration. In 
this way the amplitude may grow with energy as fast as g'^s{g'^ In sY, and there are diagrams 
doing so. In the notation of the Introduction, this means that all amplitudes are bounded 
by ~as. 

Multi-loop Feynman diagrams can usually be computed only in LLA. When they are 
summed, their leading power of In s unfortunately cancels in most t-channel colour config- 
urations. See e.g. Ref. [|T^] for concrete examples. Sometimes many subleading powers 



are cancelled as well. In the case of QED with multi-photon exchanges, all Ins powers 
are cancelled. This simply means that a non-zero sum can be obtained for these colour 
configurations only when individual Feynman diagrams are computed to the appropriate 
subleading-log accuracy, which is generally quite an impossible thing to do. 



This disaster is avoided in cut diagrams [0. This is so because all cancellations to 



occur have already taken place in building the individual cut diagrams. At a cut line, the 
Feynman propagator (|) is replaced by the Cutkosky propagator —27!-i6{K^). If a In s factor 
is to occur in the loop involving this Feynman propagator, it comes from its singularity at 
K^ = through the integral [^ 

/ ^^^_^ln(s/A2), (7) 

where A^ is determined by a mixture of the other scales (rra^, A^, Q"^) in the problem, and can 
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be either positive or negative. In LLA the In A^ factor is dropped so its precise dependence 
becomes immaterial. When this Feynman propagator is replaced by a Cutkosky propagator, 
the integral becomes a constant, and the Ins factor disappears. This corresponds to the 
cancellation of In s factors when Feynman diagrams are summed, but here the cancellation 
has already taken place once the Feynman propagator is changed into the Cutkosky propa- 
gator, even before the high-energy limit is calculated. This is why it is sufficient to calculate 
each cut amplitude in LLA. In contrast, if we first compute the high-energy limit from in- 
dividual Feynman diagrams before summing them, then cancellation occurs afterwards and 
each diagram has to be computed to a subleading-log accuracy. 

By symmetry the same thing happens on the lower parton line. Since a Ins factor is lost 
for each cut on the upper parton line, or the lower (but not both), the g and s dependence 
of a nonabelian cut amplitude with i loops, mi — 1 cuts on top, and 1712 — 1 cuts at the 
bottom, is bounded by 

/^s(/lns)^-^+\ (8) 

where M = raax(mi,m2). In other words, it is bounded by ~a^. Since 7711(7712) is also 
the number of adjoint-colour Reggeons emitted from the upper (lower) parton (see the 
discussion at the end of Sec. 3), bounds with different M apply to different Reggeon (or 
colour) channels. It is this ability of the cut diagram to extract small contributions from 
large-colour channels that makes it so valuable for unitarization. With this extraction we 
can now afford to throw away even smaller contributions at the same colour. This is what 
we mean by the extended leading-log approximation (eLLA). It differs from LLA in that it 
keeps the leading contribution at every colour, no matter how small that is. 

With eLLA and t-channel colour conservation, Reggeon number is conserved. In other 
words, we can ignore diagrams with mi 7^ 7712- This is so because the colour being exchanged 
must be contained both in the product of m,i adjoint colours, and the product of m,2 adjoint 
colours. If m,i > m,2, say, then there is another diagram, with m,[ = m2, that possesses 
all these colours that are allowed, but with an amplitude ~Q;™^ which dominates over the 
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present one whose amplitude is ^a^ ~a™^. Hence diagrams with mi 7^ m2 can be dropped 
in eLLA. 

This cut amphtude, with magnitude ~a*^ and generated by the exchange of M adjoint- 
colour objects, is just the M-Reggeon (fragment) amplitude. When all the fragments are 
summed up it becomes the usual M-Reggeon amplitude. In particular, this shows that the 
most dominant comes from a one-Reggeon exchange and is ~ a^ , agreeing with what we 
know by direct calculation |^. 

We shall discuss in the next section the factorization of Reggeon (fragment) amplitudes. 

V. FACTORIZATION, EXPONENTIATION, AND PHASE SHIFTS 

We shall consider in this section factorization of Reggeon (fragment) amplitudes into 
irreducible amplitudes. This factorization enables all amplitudes to be summed up into an 
eikonal form, and phase shift expressed as the sum of all irreducible Reggeon (fragment) 
amplitudes. 

We will call a Feynman diagram reducible, if it falls into k > 1 disconnected parts after 
the two energetic-parton lines are removed. For example, both Fig. 3 and Fig. 4 are both 
reducible. They have k = 3 irreducible parts if each of a,b,c is irreducible. A similar 
definition will be used for cut diagrams, except that all consecutive uncut gluon lines must 
be merged into single Reggeon lines before the two energetic partons are removed. For 
example. Fig. 5 has three irreducible parts but Fig. 6 has only two. 

Every member in a permuted set of Feynman diagrams has the same number k of irre- 
ducible components, but this is not so in a permuted set of cut diagrams as evident from 
the example above. 

In every permuted set there is at least one member diagram which is uncrossed. This 
means the gluon lines from different irreducible components do not cross one another. For 
example. Figs. 3 and 5 are uncrossed diagrams, and Figs. 4 and 6 are crossed. For reducible 
diagrams with k distinct components, there are actually k\ uncrossed diagrams obtained by 
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permuting the relative positions of these components. 

To discuss factorization we must first discuss how to construct and to organize the cut 
diagrams in a systematic way. We will do so by induction on the order n, assuming that we 
already know how to do it for all orders less than n. 

In this connection the remark made in the last paragraph of Sec. 3 about the labelling 
of gluon lines is relevant. The labelling can be assigned arbitrarily for one diagram in a 
permuted set, after which the location of cuts in any other diagram of the set is completely 
determined. We shall choose our labellings to facilitate factorization. 

Consider any permuted set of diagrams, li k = 1, there is no question of factorization 
so we just label the gluon lines anyway we want to. If A; > 1, we will choose to fix the 
labelling from an uncrossed diagram, such as Figs. 3. If there are n^ gluons in part a, n^ 
in part b, etc., then label the gluons in part a by the numbers 1 to n^, and those in part 
b by the numbers Ua + 1 to Ua + rib, etc. This labelling ensures cuts to occur between any 
two irreducible components, as in Fig. 5. The labelling within each irreducible component 
is assumed to have been fixed before, and it is correlated with the location of cuts inside 
these components. Fig. 5 shows one labelling that gives rise to the cuts shown. 

With cuts between irreducible components, the amplitude of the uncrossed cut diagram 
factorizes in the impact-parameter space, 

Ais,b) = -2isY[[ihj{s,b)Gj] , (9) 

j 

where Gj is the colour factor of the jth irreducible component. The order of Gj in the 

product follows the order of the irreducible components in the uncrossed diagram. To see 

why we have to go to the impact-parameter space for factorization, and to understand the 

detailed numerical factors appearing in (|^), let us briefly review the Feynman rules used to 

construct the T-matrix. Feynman propagators are taken to be M /{—q^ + m? — ie), with 

A/" = 1,771 + 7 ■ g, —g^"^ respectively for scalars, spinors, and gauge bosons. Vertex factors 

are taken directly from the coefficients of the interacting Lagrangians. Then all the i's and 

(27r)'s of the T-matrix are contained in the factor [—i / {2tt)'^Y , where i is the number of 
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loops in the diagram. In other words, we have a factor —2/(27?)^ per loop. Let q be the 
loop momentum of any loop containing a cut on top and a cut at the bottom. The loop 
integration for the impact-parameter amplitude, taken the top and the bottom Cutkosky 
propagators into account, becomes 

The factor 2s in the denominator of (|10|) is cancelled by the factor 2s from a pair of vertices. 



Since there is one less loop than pairs of vertices present, an extra factor 2s remains, which 
is shown in (P) in front of the product sign. This leaves the factor i per loop in (|10|), 



which accounts for all the i factors in (P). Finally, transverse momentum conservation links 
the irreducible components together, preventing factorization to occur. But if we Fourier- 
transform it into impact-parameter space, then the exponential factors exp(iqj ■ b) separate, 
enabling factorization to take place as shown in (|^). 

Once the labelling of gluon lines is fixed this way in the uncrossed diagram, the location 
of cuts in any other cut diagram within the permuted set is completely determined. A 
diagram in this set may come from permutation of lines within irreducible components, in 
which case it involves only lower-order diagrams and by the induction hypothesis we do not 
have to worry about it. That leaves diagrams whose gluon lines from different components 
cross one another, as in Fig. 6. Such crossed lines always fuse together some irreducible 
components and reduce their total number from A; to a number k' < k. If A;' = 1, this is 
simply a new irreducible component. If A;' > 1, then each irreducible component is of order 
< n, so it must have been included already in a lower-order consideration. In any case, 
the factorization formula (j^) is still valid, so every cut amplitude can be factorized into a 
product of irreducible Reggeon amplitudes. 

A permuted set of Feynman diagrams with k irreducible components has k\ uncross 
diagrams, corresponding to the different ordering of the components. We could have fixed 
the labelling using any of these. For example, we could fix the labelling from Fig. 7(a), 
getting the cuts shown there, (the thick lines are Reggeon hues, i.e., a group of gluon lines 



with no cuts between them), or we could have fixed the labelhng from Fig. 7(b). The 
former factorize into something proportional to GfeGaGc, and the latter is proportional to 
GaGbGc- Since the colour matrices do not commute these two are not the same. Which one 
should we use? It turns out that within eLLA it does not matter. Either will do and they 
give effectively the same amplitude (^), because their difference is an amplitude that can 
be neglected in eLLA. This is so because the difference is proportional to the commutator 
[Ga,Gfe]. Using Jacobi identity, commutator of two nested multiple commutators can be 
written as sums of nested multiple commutators. So the difference of Figs. 7(a) and 7(b) is 
given by cut diagrams with one less Reggeon line, say M — 1 instead of M, which means 
that the spacetime amphtude ~a*^ of these diagrams is negigible compared to the leading 
contribution r^af^~^ for amplitudes with M — 1 reggeons. 

Summing up all permuted sets of cut diagrams of all orders, we obtain the complete 
impact-space amphtude in eLLA to be 



A(s, h) = 2is 



1 — exp ^ ihj{s, b)Gj 



(11) 



where the sum is taken over all irreducible Reggeon fragment amplitudes. By comparing 
this with the formula in (|I]), we deduce immediately the formula for phase shift to be 

2S{s,b) = j:h,{s,b)G,. (12) 

We may combine all the irreducible amplitudes with the same t-channel colour a to write 
the phase shift as 

26{s,b) = Y.da{s,b)Ca, (13) 

a 

where Cq is the colour factor for irreducible colour a. For example, Cg = t^ x t!^ is the 
adjoint-colour factor, where t'^ and t'^ are the colour matrices of the upper and the lower 
particles respectively. Any diagram of the type Fig. 8(a), for example, will have a colour 
factor proportional to Cg of Fig. 8(b). We shall also use Ci = 1 x 1 to denote the singlet 
colour factor. Let Ma be the smallest integer required for irreducible colour a to appear in 
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the product of Ma adjoint representations (Mg = 1, Mi = 2, etc.). Then da{s, h) ~a^'", and 
contribution to this colour from larger number of Reggeons can be omitted. 

For quark-quark scattering only adjoint and singlet colours can be exchanged, its phase 
shift is therefore given by 

5(s, b) = 6s{s, b) + 5i(s, b) = ds{s, 6)t; X t;; + di{s, 6)1 X 1. (14) 

The scattering amplitude can now be computed from (|l]), and the total cross section from 
(I). To obtain the latter it is necessary to project out the colour-singlet amplitude Ai{s,0) 
of A{s,0). The necessary algebra is carried out in Appendix A, with an answer given by 
eq. ([A13| ). From that one can obtain the total cross section to be 



aris) = 2jd^hll - ^ 2e-2K+4/3) cos (2rff + 2</3 

+e-2(<-24/3)cos(2df-</3) 



(15) 



where df and df are respectively the real and imaginary parts of the phase shift di{s,b). 
For gluon-gluon scattering, more colours can be exchanged and more phase shifts have to 
be kept. 

We shall discuss in the next two sections more detailed expressions for the phase shifts 
in pQCD. 

VI. UNITARIZATION OF THE BFKL POMERON 

Let Asisjb) and Ai{s,b) be the adjoint and the singlet components of the impact- 
parameter amplitude A{s,b). Their respective leading contributions will be denoted by 
^g(s,6) and A'i{s,b). These can be obtained by substituting ( |14|) into (0) and expanding 
the exponential. In this way we obtain 

A'{s, b) = ^^(s, b)i'^ X < + A!i{s, 6)1 X 1, 

A!^{s,b) = -2idsis,b), 

A[{s, b) = -1i [di{s, b) + z64(s, bf) , (16) 
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where ^2 is the amount of singlet contained in a pair of adjoints: 

o;xtx=6ixi + ---. (17) 

For quark-quark scattering, ,^2 is given in (|A10| ) to be 

6 = u'u"{u' + u") = ^ (18) 

for SU{3c). 

Conversely, we can solve for the phase shifts from the leading amplitudes to obtain 

Ms,b) = -A's{s,b), 

d,{s,b)='-A[{s,b) + '-^,A',{s,b)\ (19) 

The impact-parameter amplitude A{s, b) is related to the momentum-space amplitude 
by a Fourier transform. According to ([1|), we have A{s,b) = {A{s,A))/2is, where 

{F)^j^Jd'Ae~^^'>F{A). (20) 

It is well known that the leading adjoint amplitude Ag{s, A) is given by the exchange of 
a Reggeon 0, with 

«(A) = ^NA'hiA), (21) 



where ||T0| 



The parameter /i is an infrared cutoff put in by hand. From (^) we can compute the adjoint 
phase shift to be 

ds{s,b) = '-A',{s,b) = ^(A;(., A)) = -^(A-V"(^)). (23) 

Note that 
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(4) = (A)". (24) 

Similarly, the leading singlet amplitude A[{s,A) is given by the exchange of a BFKL 
Pomeron P,^, p!0| , p!9[| . The details of the Pomeron amplitude is much less well known, even 
within LLA. For example, it is known that at extremely high energies, the amplitude 
Ai{s, A = 0) has an energy dependence of 

but its complete energy dependence at lower energies is complicated, even at A = and 
within LLA. 



VII. THREE-LOOP QUARK-QUARK AMPLITUDE 

As discussed in the last section, phase shifts can be computed from ( p!4D and (|1^). Un- 
fortunately the s and b dependences of the BFKL-Pomeron amplitude are not sufficiently 
well known to allow us to make a reliable calculation that way. To get an idea how unita- 
rization affects the energy dependence of the cross section, and to illustrate the formalism 
with concrete formulas, we discuss in this section the computation of quark-quark scattering 
phase shift to three-loop order. 

The quark-quark amplitude up to two loops can be found in the book of Cheng and Wu 
To| , where references to the original literature are given. The three-loop amplitude can be 



found in Ref. ^. 

As discussed in (P) and ([121) , phase shifts can be extracted from perturbative amplitudes 
of irreducible colour diagrams. To three-loop order, one obtains in this way from Table II 



of Ref. 19 that 



J2 hjis, b)Gj = /iiGi + /laiGai + /i22G22 = 4tL x t;;; + dil x 1, 






2 r 1 1 



{h)+v{h)--v^A''ll) + -v^A^ll) 



22 



h 
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V ( (/3) - li^'l!) ) + v' {{A'hh) - {h){A'l',) - {h 



V = — NAns, 
47r 



(26) 



where Gi, G21, G22 are the colour factors for Fig. 8(b), (c), (d) respectively. Strictly speak- 
ing, the amphtudes hj and the colour factors Gj are not the ones appearing in (p!3D, because 
all colour structures of the form Fig. 8(a) have already been turned into Fig. 8(b) in the 
definition of hj here. The function /n(A) is defined in ( p2D and the brackets (F) in (pO[). 
The last term in /i2i is of order g^^ ||19|, but it must be included to keep the exponents in 



(|15|) negative and the integral convergent. 

Gi carries one Reggeon fragment in the t-channel and G12, G22 each carries two. We 
therefore expect hi to be ^ag and /i2i;^22 to be ~ag. This is precisely what is shown in 



(Pq). Table II of Ref. |T9| also contains amplitudes for the reducible colour factors G" for 
n = 2,3,4, and G1G21. It can be verified there that their corresponding amplitudes can be 
obtained by the expansion of the exponential containing the phase shift in (|l|). 

Using (|A17|) to obtain the adjoint and the singlet projections of G21 and G22 for SU{3c), 
we obtain the phase shifts to be 



1 3 

dsi-s, b) = hi- -h2i - -h22, 



di{s,b) 



hoi + 2/i 



22- 



(27) 



The contribution to the adjoint-colour amplitude can be computed from (|l|), (|T^), (26) 



and (|27|) to be 



^ [1 - «(A) In s + i (a(A) In s)^ - 1 (a(A) In s)^' 
2sg^ 



,-"(A)] 



A2 



(28) 
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This is just the first four terms of (|2l|) when e""*^^)^"^^ is expanded. So to order g^ we see 
directly in this way that the adjoint amphtude is dominated by the 1-Reggeon exchange. 
The leading contribution to the singlet amplitude can be computed similarly. It is 



A[is, A) = ^\v (4/3 - 2A'I^) + v' (-AA'hh + 2A^/|) + ^v^h * (A'l',) 



g ,^— o zy ■ - V ^ -^ ' ^J ' 3 



1 

A2 



-a{A) 



4 



1 

A2 



s 



-a(A) 



(29) 



where 



Fi(A) * F,{A) = ^ / d'A'F^iA - A')F2(A'), (30) 

and an infrared cutoff A^ ^ A^ + /i^ should be introduced by hand as in (|2^) and (|2^) to 
simulate a hadronic size. This is the BFKL-Pomeron amplitude accurate to order g^. 

Since gi = hi,g2i = —ih2i, and (722 = —^^^22 are real, the appropriate combination used 
to compute the total cross section from (|15D becomes 

2 17 

2d{ + -di = -g2i + -^22, 

4 5 

2d[ - -di = g2i + -^22, 

2df+'^di = ^gu 

2df-^< = -^(7i. (31) 



VIII. 7*P AND OTHER TOTAL CROSS SECTIONS 

7*P total cross section will be discussed in this section. We will concentrate on its 
energy variation which is governed by the (unitarized) Pomeron and is thus universal. The 
magnitude of the cross section depends on the hadronic size so there is no way we can say 
much about it in pQCD without a model of the hadron. 

The experimental situation is as follows IHJ^. The variation with energy is consistent 
with a power growth, a^^s^Q) ~ s"''^'^\ with a Q-dependent exponent a{Q)-, which is an 
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increasing function of virtuality Q^ . Its value at Q = is consistent with the universal 



exponent 0.08 observed in all hadronic total cross sections pO |. 

Compared to hadronic cross sections, the new and interesting feature is the dependence of 
the exponent on Q. We suggest that this may be taken as evidence that unitarity correction is 
already at work at HERA energies, that this dependence on Q is a reflection of the slowdown 
of total cross section growth, needed to satisfy the Froissart bound In s asymptotically. The 
fact that the observed exponent is smaller than that calculated from the LLA BFKL Pomeron 
lends further support to this suggestion. 

We shall now explain why the Q-dependence of the exponent can be taken as sign of a 
slowing growth. To saturate the Froissart bound, the total cross section asymptotically will 
have the form aT{s,Q) = ao{Q)lD. (s/A^), where A is the scale parameter to measure the 
energy ^/s with. If Q is much larger than the other dimensional variables in the problem 
(masses, Aqcd, and A), we expect A to be determined by Q. On dimensional grounds 
the simplest dependence would be A = cQ for some dimensionless constant c. In any case 
A{Q) is expected to increase with Q. Extrapolating back to HERA energies, the cross 
section becomes aT{s,Q) = aQ{Q)f\s/ Is?{Q)), for some function f\s') which approaches 
In^ s' asymptotically. For a sufficiently small range of s, the energy variation can always be 
simulated by a power growth, with the effective exponent a = d\nf{s')/d\n{s') given by the 
slope of the curve in a log- log plot. Since this slope is positive at HERA energies and it must 
decline to zero at asymptotic energies, it is reasonable to assume it to be a monotonically 
decreasing function of s'. In other words, the rate of growth of f{s') decreases with energy. 
Now for a given range of s, the corresponding range of s' becomes smaller for larger Q, 
thus placing the data at a region of faster growth. Hence the effective exponent a{Q) would 
be an increasing function of Q, as observed. We can either take that as a prediction for 
the observed increase of a{Q) with Q, or reverse the argument and interpret the observed 
variation as an indication for the declining slope of ln/(s') with increasing Ins', in an effort 
to comply with the Froissart bound at asymptotic energies. 

Qualitatively, this prediction of a{Q) increasing with Q is quite robust, as it is quite 
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independent of the detailed form of f{s') and A(Q). If we know the functions / and A then 
the prediction can also be verified quantitatively. Knowing A{Q) we can convert the observed 
aT{s,Q) into a function of s' and Q. Now plot In ctt against Ins' for every (sufficiently 
large) Q, and also ln/(s') against Ins'. The prediction then says that by a suitable up-down 
movement of the experimental curves, every one of them can be made to fall on the universal 
curve ln/(s') (moving an experimental curve up and down is equivalent to adjusting the 
unknown function o"o(Q)). 

Are f{s') and A{Q) calculable in eLLA? In principle / is, but unfortunately A is not. It 
is in the nature of a leading-log approximation, extended or not, that only the highest power 
term of Ins is kept. That means even if we knew A as in ln(s/A^) = Ins — InA^, it would 
have been dropped during the course of the calculation. So unless we can go beyond eLLA, 
A{Q) must be taken as a parameter. As to /(s'), which is the same as /(s) in eLLA, it is 
given essentially by ([T5|) and ([191) . These formulas compute quark-quark total cross section, 
not 7*P. However, energy variation is governed by the Pomeron and should be universal, so 
up to the unknown function o"o(Q) we may simply take them to be the same and take (|T3|) 
to be /(s) for the purpose of the test above. 

Even with an unknown A{Q) the prediction and the test are still not trivial, because in 
general a function aT{s,Q) of two variables cannot be fitted by two functions ao{Q),A{Q) 
of one variables. 

As mentioned before, the BFKL- Pomeron amplitude is not sufficiently well known as yet 
for the universal curve /(s') to be calculated accurately at this moment. However, we can 
illustrate this prediction by using the three- loop result to calculate /(s') approximately ^1 



A theory of massless quark-quark scattering lacks an intrinsic distance scale. This causes 



an infrared divergence which is cut off by the parameter fi in (|22D. By making a scaling 
change b -^ jjh, we see from (|I]) and (Q) that axi-s) = fi~^aj^ (s), with cr^- (s) obtained from 
axis) by setting /i = 1. Thus /i~^ is the distance unit to measure the cross sections in. It 
simulates hadronic sizes that would come in to provide an infrared cutoff in hadronic cross 
sections. It can be considered simply as a part of the unknown function (Jo{Q). In this way 
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we can use ([T2D, (0), and ( p6D with /z = 1 to calculate o"t[7*-P] and /(s') up to an overall 
normalization a^lQ). 

For the parametric function A{Q) we shall take a simple form cQ + Aq. This is simply 
cQ for large enough Q, but since we do not know where the cutoff of Q is, we put in this 
parameter Aq to accommodate the smaller-Q data. With c = 4, Aq = 0.2, and a properly 
chosen ao{Q), the computed three-loop result is shown as a solid curve in Fig. 9. They 
agree quite well with the experimental data. In particular, the Q dependence of the energy 
growth exponent a{Q) is reproduced. The dotted curves show an energy variation of s°°*, 
appropriate for hadronic total cross sections. They are placed there to show that the rate 
of energy growth for the theoretical calculation grows with Q^, and demanded by the data. 
The dash curve is s°'^, as given by the BFKL Pomeron. 

We can use the fitted values of ctq{Q), c, and Aq to compute a{s, Q)/croiQ) as a function 
of s' = s/ h?{Q). This is shown in Fig. 10 together with the three-loop energy function f{s') 
(solid curve). 
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APPENDIX A: COLOUR DECOMPOSITION 

The necessary colour decompositions for quark-quark scattering are worked out in this 
appendix. 

The U{Nc) colour matrices ta (0 < a < N^ — 1) are conventionally normalized to be 

Tr(t,t^) = -6a/3. (Al) 

This leads to the completeness relation 
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{ta)ij{ta)kl — -^^il^kj, (A2) 

where a sum over the repeated index a is understood. We shall use the Latin indices, 
a = a {1 < a < N^ — 1), to label the SU{Nc) generators. The remaining U{Nc) generator 



is to = l/v2iVc. To simplify writing, we shall use (• ■ ■) to denote the trace Tr(- ■ ■), and the 
index a to represent the generator ta- In this notation, the structure constant ifais-y defined 
by the commutation relation 

{taitp] = ifap-yt-y (A3) 

is given in terms of the traces of the generators to be 

i/„^^ = 2(a/37-a7/5). (A4) 

It is clear that fap-^ is totally antisymmetric in its indices, and that /0/37 = 0. 
The following formulas follow from the completeness relation (|A2|) : 

{Aa){Ba) = ]^{AB), 

{Aa){Ba) = \{AB) - ^^{A){B), 

{AaBa) = ]^{A){B), 

{AaBa) = \{A){B)-^^{AB). (A5) 

They will be used to compute colour decompositions. 

Let t'^ and t'^ denote the colour matrices of the upper quark and the lower quark, respec- 
tively. A colour diagram like Fig. 8(a) has the colour factor Gi = t^ x t'!^. The colour factor 
for the reducible diagram when this is repeated n times is 

Gi = ti^4 ■ ■ ■ ti„ X titi . . . t'i ^ (E„),.t; X c (A6) 

The last expression comes about because the only colour allowed to be exchanged between 
two quarks is either a colour singlet 1x1, or an adjoint colour t!^ x t'^. The coefficients 
{En)^u can be computed from (|A1|) to be 
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{En)f,u =4(/iaia2---a„)(i/aia2---a„) = An{fiu) + Bn{fi){i^). 



(AT) 



Using (|A5|) to contract the pair of indices a„, we obtain a recursion relation for An and Bn'- 



Bn = - (a„_i - —Bn-, 



(A8) 



This pair of equations can be diagonahzed using the variables C^ = y4„ ± _B„, and solved to 
obtain 



^.--4(iT^)e.-±4 



±^ IT 



Nr 



Subsituting this into (|A6|) and ( |A7|) , we obtain 



^" = i (f^" + ^") ^^1 + ^^"^'- ^ *^ 



m'm" 



/\n— 1 



M 



//\n— 1 



-M 



ixi + [Kr-(-«T]tLx<, 






M 



iV. 



From this it follows that 



(A9) 



(AlO) 



exp i2id8t'^ X t;;;) = [^''exp {2iu'ds) + w'exp {-2iu"ds)] 1x1 
+ [exp {2m ds) - exp {-2iu"ds)] t'^ x t;;;. 



(All) 



For SU{3c), u' = 1/3 and u" = 2/3, so ( [ATI] ) becomes 

exp (2^d8t:, X t;:) = - [2e2^'^«/=^ + e-^""'/^ 



1x1 



+ 



g2id8/3 _ g-4id8/3 



tL X <. 



(A12) 



The colour-singlet component of the impact-parameter amplitude Ai{s, h) is therefore given 
by 

1 



A(s,&) = l-- 



2g2i(di+dg/3) _^ g2J(di-2a!8/3) 
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(A13) 



Separating the phase shifts into their real and imaginary parts, di = df + idl, and using 
we finally obtain the formula for total cross section in terms of phase shifts to be 



aT{s) = 2Jd''hh-^- 



+e-2(4-24/3)cos(2rff-4df/3 



(AM) 



For the purpose of calculating phase shifts from perturbation theory, other colour de- 
compositions are required, especially those with two t-gluon lines in the colour diagram. Let 
G2„ be the colour factor for a colour diagram with 2 t-gluons and n 'horizontal' gluons. For 
example, the colour factor for Fig. 8(c) is G21 and the colour factor for Fig. 8(d) is G22- 
The colour factor of a triple gluon vertex in these diagrams is if aha read counter-clockwise, 
so that 



G2I — {iface){ifedb)tatb X i^c^d) 

G22 = {ifahg){ifceh){iffed){iffbg)t'at'b X t%. 



(A15) 



We will now show that 



G2n = \K-'[Nl-l)l^l-\(^T\'^^il 



In particular, for SU{?>c), 



G 



2 1 
21 = -1x1- -ti X t'', 

3 2 ^ ^ 



(A16) 



3., 



G99 = 2ixi--tl X r. 



(A17) 



Eq. ( |A16|) can again be obtained by induction, as follows. First of all, define In{oi[3\'^5) 



by 



G2n = In{aP\l6)t% X t'^^tl 



(A18) 



where /„ contains n pairs of triple-gluon colour factors [if...)- See Figs. 6(c),(d),(e). Strictly 
speaking, the subscripts oft' and t" should be Latin indices and not Greek, but since /q.. = 0, 



we may extend the subscripts to Greek indices as shown in (|A18|) for n > 1. For n = 0, we 
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should identify G20 with Gi. Next, as only singlet and adjoint colours can be exchanged, 
we can write 

G2„ = {F„)^,t'^ X t'l = Pnlxl + Q^t', X t^, 
(F„)^, = 4{fia(3)In{a(3\^6){u^6). (A19) 

Now write a recursion formula for In{al3\'~f6){i"~f6), as indicated in Fig. 8(f): 

Iniap\'jS){u'jS) = In-li(yP\pCr)iif-frp){ifar5){^lS) 

= In-i{aP\pa){u[T,p][a,T]). (A20) 

This is valid for all n > 1, though we should note that G20 = G^, hence /o(a/?|pcr) vanishes 
unless all four are SU{Nc) indices, in which case Io{ab\rs) = dar^bs- The repeated r index 
can be summed up using (|A5|) to get 



(z.[r,p][a,r]) = \ {-{p){ua) - {a){up) + {u){pa) + {l){upa)) . (A21) 

The first two terms do not vanish only for z/ = and a = 0, but then since /q.. = 0, these 
two terms never contribute. This leaves the last two terms. It is now convenient to consider 
separately z/ = and u = n. In the first case, the last two terms are identical and the right 



hand side of (|A21|) becomes JNc/2{pa) = JNc/86p„. Now In~i{al3\pa) contains {n — 1) 



pairs of triple-gluon colour factors. Its bottom pair takes on the form {ifp(t,(;){if^^a)- When 
contracted with 6pa, we get {ifp^(^){ifi^,f,p) = A^c^c^- Hence 

In-i{a(3\pa)6p^ = Ar,/„_2(«/?|CO'5c« = " " " = K"'lo{al3\pa)6p^. (A22) 

Substituting this back into eqs. ( |A2CI|) and (|A19| ), one gets 



(F„)^o = Mpa(3)./Nj8Nr'lo{a(3\pa)6p 



rn+i I I 1 



This expression of P„ agrees with the one quoted in ( [A16| ) 
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^ ^00 = ]n: (l-^). (A23) 



Now we return to ( |A21|) and consider the case when v = n. In that case the right hand 
side of that equation becomes Nc{npa), so we obtain a recursion relation from ( [A19|) to be 

1 



y^njiin „ -' 'cl^-f^n— 1 j/^n- 



(A24) 



Since the adjoint component of G20 = G^ is given in ( JAIOD to be Qq 
—1/Nc, we can solve (|A24|) to obtain 



{u'Y - {-u") 



li\2 



Qn = T^NcQn 



N, 



Qi 



1 rNr 



2 V 2 



n-1 



(A25) 



agreeing with the expression given before. This completes the proof of the formula quoted 



before in (|A16|) . 
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FIGURES 
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FIG. 1. A tree diagram depicting n quanta of momenta ki emitted from an energetic particle. 
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FIG. 2. A high-energy scattering diagram with cm energy -^s and momentum transfer A. 
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FIG. 3. An uncrossed Feynman diagram with three irreducible components. 




FIG. 4. A crossed Feynman diagram with three irreducible components. 
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FIG. 5. An uncrossed cut diagram with three irreducible components. 




FIG. 6. A crossed cut diagram with two irreducible components. 
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FIG. 7. Two uncrossed cut diagrams differing from each other by the permutation of the first 
two components. The thick lines are Reggeon fragments. 
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FIG. 8. Various colour diagrams referred to in the text and in Appendix A. 
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FIG. 9. 7*P total cross section as a function of photon virtuality Q^ and cm. energy y/s. 
The sohd curve is a fit of the three-loop eq. (|26|), with h.{Q) = [0.2 + 4Q]GeV and fJ.{Q) chosen 
arbitrarily. The dotted and the dashed curves represent an s^'^^ and an s^'^ variation respectively. 
Data are taken from Ref. [01. 
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FIG. 10. The universal energy variation curve /(s') obtained from the three-loop formula (| 
with arbitrary normalizations. The data are those of Fig. 6, with the same A((5), but with all low 
energy points with s < lOOO(Gey)^ removed. 
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